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ABSTRACT 

We present an analysis of the far-infrared (FIR) spectral energy distributions (SEDs) of two massive 
K-selected galaxies at z = 2.122 and z — 2.024 detected at 24/im, 70^m, 160/im by Spitzer, 250/ini, 
350/im, 500/ini by BLAST, and 870/ini by APEX. The large wavelength range of these observations 
and the availability of spectroscopic redshifts allow us to unambiguously identify the peak of the 
redshifted thermal emission from dust at ^ 300/im. The SEDs of both galaxies are reasonably well fit 
by synthetic templates of local galaxies with L/jj ~ lO^^L© - IQ'^^Ijq yet both galaxies have Ijir ^ 
1O^^L0. This suggests that these galaxies are not high redshift analogues of the Hyper-LIRGs/ULIRGs 
used in local templates, but are instead "scaled up" versions of local ULIRGs/LIRGs. Several lines 
of evidence point to both galaxies hosting an AGN; however, the relatively cool best fit templates 
and the optical emission line ratios suggest the AGN is not the dominant source heating the dust. 
For both galaxies the star formation rate determined from the best-fit FIR SEDs (SFR(L/fl)) agrees 
with the SFR determined from the dust corrected Ha luminosity (SFR(Ha)) to within a factor of ~ 
2; however, when the SFR of these galaxies is estimated using only the observed 24/im flux and the 
standard luminosity-dependent template method (SFR(24/im)), it systematically overestimates the 
SFR by as much as a factor of 6. A larger sample of 24 K-selected galaxies at z ~ 2.3 drawn from the 
Kriek et al. (2008) GNIRS sample shows the same trend between SFR(24/^m) and SFR(Hq!). Using 
that sample we show that SFR(24/^m) and SFR(IIa) are in better agreement when SFR(24/im) is 
estimated using the log average of local templates rather than selecting a single luminosity-dependent 
template, because this incorporates lower luminosity templates. The better agreement between SFRs 
from lower luminosity templates suggests that the FIR SEDs of the BLAST-detected galaxies may 
be typical for massive star forming galaxies at z ^ 2, and that the majority are scaled up versions of 
lower luminosity local galaxies. 



Subject headings: infrared: galaxies 



1. INTRODUCTION 

It is well known that in the local universe the obscura- 
tion towards star forming regions is correlated with star 
formation rate (SFR; e.g., Wang & Heckman 1996; Hop- 
kins et al. 2001; Buat et al. 2005, 2007), with the most 
active galaxies frequently being the most dust obscured. 
In the most extreme galaxies, the star forming regions 
can become optically-thick (e.g., Genzel et al. 1998), 
which makes mid- and far-infrared (MIR, FIR) data crit- 
ical for a complete metric of their SFR. At z ~ 2, where 
the average massive galaxy is forming stars at rates ~ 
100 - 1000 times higher than in the local universe (e.g., 
Perez-Gonzalez et al. 2005; Juneau et al. 2005; Damen 
et al. 2009), MIR- and FIR-determined SFRs are critical 
for understanding the buildup of today's massive galax- 
ies. 

The advent of the MIPS instrument onboard Spitzer 
has facilitated the first deep photometry at 24/im and 
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numerous studies have used these data to infer the star 
formation history of the universe up to z 2 (e.g., Le 
Floc'h et al. 2005; Perez-Gonzalez et al. 2005; Caputi et 
al. 2007; Reddy et al. 2008; Magnelli et al. 2009). These 
studies suggest that obscured star formation dominates 
in massive galaxies at z 2; however, the conversion 
from observed 24/im fiux to a total infrared luminosity 
(L/fl), and subsequently a SFR at z ~ 2 requires signifi- 
cant extrapolation based on local FIR templates. 

The most common method for determining the L/jj 
of a galaxy from 24/im photometry is to artificially red- 
shift local luminosity-dependent MIR/FIR templates to 
the redshift of the galaxy and then choose the L/^j of 
the template that predicts a 24/im flux closest to the ob- 
served 24/im flux. By applying this method, or more so- 
phisticated but related versions of it, some authors have 
found that the local templates produce estimates of L/j? 
and SFR(L/fl) for distant galaxies that are in agreement 
with other SFR indicators (e.g., Elbaz et al. 2002, Mar- 
cillac et al. 2006, Reddy et al. 2006; Daddi et al. 2007a; 
Papovich et al. 2009; Kartaltepe et al. 2010), while oth- 
ers have argued that this method tends to systematically 
overestimate the I^ir and hence the SFR (e.g., Papovich 
et al. 2007; Rigby et al. 2008; Murphy et al. 2009). 
Some of these differences can probably be reconciled by 
the different selection criteria (e.g., Reddy et al. 2010), 
redshift range (e.g.. Murphy et al. 2009) and luminosities 
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(e.g., Papovich et al. 2007) of the above samples; how- 
ever, the rehabihty of observed 24/im fluxes as a metric 
of the Lju and SFR over the full range of galaxy masses 
and SFRs at z ^ 2 remains unclear. 

The next step will be to directly observe the cold dust 
SED, rather than extrapolate it from the rest- frame 6 - 
10/im region, which can be complicated by a superpo- 
sition of Polycyclic Aromatic Hydrocarbon (PAH) emis- 
sion lines and silicate absorption features. The peak of 
the cold dust emission in FIR luminous galaxies in the 
local universe occurs between 40/im to 150/im which cor- 
responds to 120/xm to 450/im at z ~ 2. In the next few 
years the Herschel telescope will provide unprecedented 
amounts of data at these wavelengths; however, the re- 
cent flight of the BLAST telescope (Devhn et al. 2009) 
has already provided some of the first deep observations 
of distant galaxies at 250/im, 350/xm and 500/xm. 

In this paper we take advantage of the BLAST ob- 
servations of the Extended Chandra Deep Fields South 
field (ECDFS) to directly measure the FIR SEDs of z - 
2 galaxies. Two star forming galaxies in the Kriek et al. 
(2008) NIR spectroscopic sample are located in that field 
and are detected or have upper limits at 24/im, 70/im, 
and 160/im in the FIDEL (Magnelh et al. 2009) and 
SWIRE (Lonsdale et al. 2003) surveys, and are also 
detected by BLAST (Dye et al. 2009). The ECDFS 
was recently surveyed by the LABOCA instrument on 
the APEX telescope, and both galaxies were detected at 
870/im (Weiss et al. 2009). By combining these data 
sets we have produced some of the first FIR SEDs that 
sample the cool dust bump of z ~ 2 galaxies with good 
resolution and allow us to constrain their L//j directly 
from the FIR SED. Throughout this paper we assume a 
Ho = 70 km s"i Mpc"\ Vim = 0.3, r^A = 0.7 cosmology 
when computing luminosity distances. 

2. GALAXY SAMPLE AND DATA 

The NIR spectroscopic sample of Kriek et al. (2008) is 
a magnitude-limited (K < 19.7 Vega) sample of 36 galax- 
ies with a median z ^ 2.3. The sample is dominated by 
massive galaxies (M > lO^^M©), approximately 50% of 
which have detected emission lines. Two of the galax- 
ies in the sample are extremely bright at 24/im (S24/im 
> 0.5 mJy), and both are located in the ECDFS (ID#'s 
ECDFS-4511 and ECDFS-12514 from Kriek et al. 2008). 
Both galaxies have strong Ha emission and Kriek et al. 
(2008) used this to determine spectroscopic redshifts of 
z = 2.122 and z = 2.024. The ID numbers of the galaxies 
in the BLAST survey of the ECDFS are BLAST J033138- 
274122 and BLAST J033242-275511. They are also de- 
tected at 870/im in the 200-hr map from the LABOCA 
array on the APEX telescope (Weiss et al. 2009, ID#'s 
LESS J033139.6-274120, and LESS J033243.3-275517). 

Using the vO.5 pubhc-release FIDElJ3 MIPS imaging 
data of the ECDFS we have measured total 24/im and 
70/im fluxes for these galaxies. For the 24/(m data the 
fluxes were measured within 11" diameter apertures us- 
ing the SExtractor package (Bertin & Arnouts 1996) and 
an aperture correction from the MIPS data handboolsQ 
was applied. For the 70/im data we used the 24/im cen- 
troid and measured fluxes within 32" apertures and ap- 

http://irsa.ipac.caltech.edu/data/SPITZER/FIDEL/ 
http:/ /ssc. spitzer.caltech.edu/mips/dh/ 



plied the aperture correction from the MIPS data hand- 
book. The ECDFS was observed at 160/(m as part of the 
SWIRE survey (Lonsdale et al. 2003). From the data we 
measured 160/im fluxes within 64" apertures, again us- 
ing the 24/im centroid. The SWIRE 160/im data is much 
shallower than the FIDEL data; however, ECDFS-12514 
is detected at the 1.5a level (a similar detection was made 
byDyeetal. 2009 using the SWIRE data). ECDFS-4511 
is not detected; however, we calculate a 2(7 upper limit 
based on empty aperture fluxes from the SWIRE pho- 
tometry. 

In our SED fitting we use our measured 24/tm, 70/fm 
and 160/im fluxes, the 250/(m, 350/(m, and 500/im fluxes 
measured by Dye et al. (2009), and the observed 870/im 
fluxes from Weiss et al. (2009). For ECDFS-4511 we 
do not use the 160/im upper limit in the fitting. These 
fluxes are listed in Table 1. 

A joint analysis of the FIDEL 24/im and BLAST FIR 
data by Chary & Pope (2010) has shown that source 
blending is a problem for many of the BLAST detections. 
We have used the 24/im data to test for contamination of 
the longer wavelength data of these galaxies by blended 
sources. 

The galaxy ECDFS-12514 is well-isolated on the 24/(m 
image. The nearest 24/im source is 10" away; however, 
it is extremely faint, only 4% of the flux of the main 
galaxy. The only nearby bright 24/im source is 21" away 
which means if it was luminous at any of the longer 
wavelengths, excluding 500/im, it should be resolved as 
a separate source. That galaxy is only 15% of the main 
galaxy flux and shows no evidence of a detection in any 
of the longer wavelength channels, suggesting that the 
70^'870/im measurements for ECDFS-12514 are proba- 
bly not contaminated by blended sources. 

Galaxy ECDFS-4511 is not completely isolated. There 
are two well-detected 24/(m sources within 30" . The clos- 
est is 7" away and is 15% of the flux of the main galaxy. 
The next closest is 11" away and is 37% of the flux of 
the main galaxy. If these companions are luminous in 
the FIR, neither would be resolved in the FIR data and 
could cause the FIR flux to be overestimated. Given 
that at z ^ 2, the SED of the main galaxy has been 
redshifted to where the dust SED peaks in the BLAST 
bands and the observed 24/im/FIR flux ratio is maxi- 
mal, we expect that the worst-case scenario would be if 
the nearby galaxies were at the same redshift and had 
similar FIR SEDs as the main galaxy and therefore con- 
tribute equally to the FIR flux. If so, it suggests the FIR 
fluxes and \jir could be overestimated by as much as ~ 
30%. If the companions are at either higher or lower red- 
shift, they are likely to have higher 24/iin/FIR flux ratios 
and therefore the overestimate is likely to be less. As we 
discuss below, this possible overestimate of the L/a of 
ECDFS-4511 does not change any of the conclusions of 
this paper. 

3. TOTAL INFRARED LUMINOSITIES AND STAR 
FORMATION RATES 

3.1. Spectral Energy Distributions and Fitting 

In the top panels of Figure 1 we plot the observed U 
870/im flux densities of the galaxies. The 24/im 
870/im flux densities, which are from dust emission are 
plotted as black points and the U— >8/im flux densities, 
which are from stellar photospheres are plotted as blue 
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Fig. 1.— Top panels: U 870Atm SEDs for galaxies ECDFS-4511 (left) and ECDFS-12514 (right). The TJ^Sfim photometry is plotted 
in blue, and the best fitting BC03 model to those data is plotted in black. The FIR photometry from Spitzer, BLAST and APEX are 
plotted in black. The best fitting CEOl model to those data is shown in red. The grid for the FIR SEDs is shown in the bottom left 
panels. The dot dashed line in the upper panel shows the of the local template with the same L/^ as the best-fit template. Both 
galaxies are inconsistent with the local template and are best described as "scaled up" versions of lower luminosity galaxies. The dotted 
FIR SED in the main panel shows the CEOl template that would be chosen based on only the observed 24^m flux. For these galaxies 
those templates overestimate the L/^ by factors of 3.0 and 6.2, respectively. The bottom right panels show 3" X 3" NICMOS images of 
the galaxies. 



TABLE 1 
MIR/FIR Photometric Data 



Wavelength ECDFS-4511 ECDFS-12514 Survey/reference 

fiux (mjy) flux (mjy) 

24 fim 0.57 ± 0.02 1.42 ± 0.03 FIDEL, this paper 

70 fim 1.6 ± 1.2 4.2 ± 1.8 FIDEL, this paper 

160 iim < 40 34 ± 20 SWIRE, this paper 

250 48 ±11 35 ± 11 BLAST, Dye et al. (2009) 

350 urn 44 ± 9 44 ± 8 BLAST, Dye et al. (2009) 

500 fim 19 ± 6 24 ± 6 BLAST, Dye et al. (2009) 

870 ^tm 5.9 ± 1.3 4.7 ± 1.2 LESS, Weiss et al. (2009) 



Note. — 1) All fluxes are measured fluxes and do not include statistical deboosting corrections. 2) ECDFS-4511 has two 24/im sources 
within 11" of the main galaxy. These galaxies cannot be resolved in the FIR data. We estimate that if they have substantial FIR flux they 
could inflate the FIR fluxes, and L/j; by as much as 30% (see text). 



points. For reference, the best fit Bruzual & Chariot 
(2003, hereafter BC03) model to the U— T'S/ini data from 
Muzzin et al. (2009) has also been plotted as the solid 
black line. Figure 1 shows that there is an unambiguous 
peak in the FIR SEDs of the galaxies at ~ 300/xm, a clear 
signature of the redshifted thermal emission from dust. 

We fit the 2Afim 870/xm data using the models of 
Chary & Elbaz (2001; hereafter CEOl). The CEOl mod- 
els are a grid of 105 synthetic SEDs that were constructed 
to reproduce the correlations between the observed MIR 



and FIR fluxes of local galaxies. Each model SED is as- 
sociated with an Ljji ranging between 3 x 10* L0 and 
6x 10^^ L0. We fit the 24/im -> 870/im data using each of 
the 105 template SEDs using a minimization method 
where the amplitude of the SED is the only free param- 
eter. 

The galaxy ECDFS-4511 is the most FIR-luminous of 
the two galaxies and is well-fit by the SED template of a 
ULIRG. In the CEOl models that template is associated 
with an Lir = l.lxlO^^ l^; however, ECDFS-4511 has 
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L/fl = 8.9^]^'4X 10^^ L0, eight times larger than the lu- 
minosity of the local template. 

The second galaxy, ECDFS-12514, is three times 
brighter at 24/im than ECDFS-4511; however, it has a 
lower Ljji. The best- fit template is a LIRG template 
with strong PAH features and is associated with Lju = 
1.2x10" Lq. The LiR of ECDFS-12514 is 5.7lJ;iXl0i2 
Lq, almost 50 times more luminous than the local tem- 
plate. If the local template is correct, the large 24/im flux 
from ECDFS-12514 likely arises from strong rest-frame 
7.8/im PAH emission redshifted into the observed 24/im 
bandpass. 

The key result in Figure 1 is that although the z — 2 
galaxies have FIR SEDs that are similar to local FIR lu- 
minous galaxies, their luminosities are roughly an order 
of magnitude larger. It appears that both galaxies are 
strongly "scaled up" versions of much lower luminosity 
galaxies in the local universe. The dot-dashed line on 
the grids in Figure 1 show the of the fit to the 
local template that is associated with the best fit Lm 
of ECDFS-4511 and ECDFS-12514. Those x^s are sig- 
nificantly worse than the best-fit and show that the 
SEDs of the z ^ 2 galaxies are inconsistent with having 
the same SED shape of local templates of similar lumi- 
nosity at > 25(7. 

3.2. Comparison of FIR and 24^m Star Formation 

Rates 

Direct constraints on the cool dust SED allow us to 
address the issue of how well L/j^, and SFR(L/fl) can be 
determined for 2: ~ 2 galaxies using 24/im data alone. 
The most common practice for determining L/^j from 
24/im photometry is to select the local template with an 
L/fl, that would match the observed 24/xm flux if red- 
shifted to the redshift of the source galaxy. In Figure 1 
we plot the CEOl template that would be selected based 
on using that method as the red dotted line. 

The extraordinary 24/im flux of both galaxies requires 
that they use the most luminous templates in the CEOl 
library. The implied \.ir for ECDFS-4511 and ECDFS- 
12514 is 2.8x10" L© and 3.8x10" Lq, respectively. 
These templates overestimate the Ijir from the full SED 
fit by factors of 3.0 and 6.2 respectively, and it is clear 
from Figure 1 that they cannot simultaneously fit both 
the 24/im flux densities and the FIR flux densities. If 
these galaxies are representative of most FIR-luminous 
galaxies in the distant universe, it suggests we may be 
systematically overestimating total SFRs by similar fac- 
tors. 

Similar results have also been found by Papovich et al. 
(2007) and Murphy et al. (2009). Papovich et al. (2007) 
compared L/jj determined from 24/7,m data to that de- 
termined using the combined 24/im, 70/im, and 160/im 
photometry and found that for the most luminous galax- 
ies (S24 > 250/iJy) the 24/im photometry overestimated 
the l^iR by factors of 2 - 10. Using a sample of 14 MIR 
luminous galaxies at 1.4 < z < 2.6 with Spitzer IRS 
spectra and SCUBA 850/im photometry Murphy et al. 
(2009) found that 24/im fluxes combined with local tem- 
plates systematically overestimated the Ljr by a factor 
of ~ 5. 

3.3. Comparison of FIR and Ha Star Formation Rates 




12 3 

A, (Stors) 



Fig. 2. — A„ determined from the Balmer decrement vs. 
determined from a fit to the galaxy stellar population for the 12 
emission line galaxies from Kriek et al. (2007), Kriek et al. (2009a), 
and van Dokkum et al. (2005) with both H/3 and Ha within the 
NIR window. Only two galaxies have secure H/3 detections, the 
remainder are plotted as la lower limits. The two BLAST-detected 
galaxies are indicated as the red lower limits. The dotted line 
shows the one-to-one relation, and the solid line is the relation E(B- 
V)ateHar = 0.44E(B-V)„e(,uiar Suggested by Calzetti et al. (1997). 
The two galaxies with H/3 detections are in reasonable agreement 
with the Calzetti relation. All of the galaxies with lower limits 
are consistent with the Calzetti relation but only about half are 
consistent with the one-to-one relation. 

The FIR-derived SFRs can also be compared to the 
dust corrected Ha SFRs for these galaxies. Ha equiva- 
lent widths (EWs) were measured for both galaxies us- 
ing moderate-resolution NIR spectra from SINFONI by 
Kriek et al. (2007). We convert their EW(Ha) to an Ha 
line flux by multiplying by the continuum flux at 6563A 
estimated using the BC03 SED model fits from Muzzin 
et al. (2009). 

We assess the extinction for the Ha flux using both 
the fit to the stellar SED (U^8/tm) and the Balmer 
decrement. For the SED fits we use the Muzzin et al. 
(2009) best-fit value of A^, which is derived from fitting 
I5C03 models using solar metallicity and the Calzetti et 
al. (2000) dust law (hereafter Ay^sEo)- For the A„ from 
the Balmer decrement (hereafter A^^Baimer) we assume 
the Calzetti et al. (2000) dust law and an intrinsic ratio 
of Ha/H/3 = 2.86. Ahhough both galaxies have high S/N 
Ha detections in the SINFONI spectra, neither has a se- 
cure H/3 detection. Therefore instead of a true Balmer 
decrement we calculate la lower limitsQ using the la up- 
per limits on H/3 determined by Kriek et al. (2007). The 
H/3 upper limits have been corrected for absorption by 
the stellar population using a BC03 SED fit to the broad- 
band photometry and the NIR spectroscopy. 

In Figure 2 we plot A^^Baimer vs. the Ay^sED for the 
two galaxies as the red lower limits. For comparison, 
measurements of the 9 other z ^ 2 galaxies in Kriek et 
al. (2007), and the galaxy from Kriek et al. (2009b) are 
also plotted. Two of these galaxies (ID#'s 1030-807 and 
CDFS-6202) have secure H/3 detections, the remainder 
are plotted as la lower limits. The one-to-one relation is 
shown as the dotted line. 

* Upper limits on H/3 translate to lower limits on Ha/H/3 and 
lower limits on A^^Balmer 
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The two galaxies with H/3 detections do not agree weh 
with the one-to-one relation. The lower limits for approx- 
imately half the galaxies are also inconsistent with the 
one-to-one relation. A similar disagreement amongst A^, 
measurements was found for local galaxies by Calzetti 

(1997) . She showed that this difference can be under- 
stood as greater extinction towards the star forming re- 
gions, where the Ha originates, than to the overall stellar 
population. The Calzetti (1997) relation for local galax- 
ies is E{B-Y),tars = (0.44 ± 0.03)E(B-V)„e6«iar. In Fig- 
ure 2 we plot that relation converted to A„ assuming 
Rv = 4.05 (Calzetti et al. 2000) as the solid line. The 
Calzetti (1997) relation seems to provide a better de- 
scription of the overall population of z ~ 2 star forming 
galaxies, albeit that most galaxies only have lower limits. 

Given the uncertainty in the Ay_Baimer, and the well- 
determined Ay^sED we have chosen to dust correct the 
Ha flux using the Ay,sED but using the E(B-V) stars — 
(0.44 ± 0.03)E(B-V)„e6«iar relation. We convert this to 
a SFR using the Kennicutt (1998) relation, SFR(H„) = 
0.57 ■ L(Ha)/l. 26x10"'*^ ergs s~^, where the factor of 
0.57 has been added to convert from a Salpeter IMF 
to a Chabrier IMF. Applying these conversions we find 
that ECDFS-4511 and ECDFS-12514 have SFR(Ha) = 
620tl% Mq yr-i and 12101^^° Mq yi-\ respectively. 
We note that the majority of the error budget for these 
measurements comes from the uncertainty in the extinc- 
tion correction, not from random errors in the measure- 
ment of the Ha line flux. 

We convert the Ljj^ to a SFR using the Kennicutt 

(1998) relation, SFR(L/fl,) = 0.57 • L/i^/5.8x IO^Lq, 
where again the factor of 0.57 is to convert to a Chabrier 
IMF. We derive SFR(L/ij) = SlOtHo Y^"^ and 
560t}g° Mq yr^i for ECDFS-4511 and ECDFS-12514, 
respectively. 

Comparing the two SFR indicators shows that they 
agree to within roughly a factor of 2 for both galaxies. 
Such an agreement is good considering that there are 
still systematic uncertainties in the conversion of both 
Ljji and L(Ha) to a SFR that are of that order (Kenni- 
cutt et al. 2009). 

It is worth noting that the agreement in SFRs only oc- 
curs because L/^ is determined from the full FIR SED 
fits, and because the extra extinction correction has been 
applied to the L(Ha). If we do not use the Lm from the 
FIR SED fit and instead use the SFR(24/im) we find 
that the IR SFR is larger than the SFR(Ha) by factors 
of 4.8 and 3.4 for ECDFS-4511 and ECDFS-12514, re- 
spectively. Likewise, if we extinction correct the L(Ha) 
using the A^^sed without extra extinction we find that 
the SFR(L7;j) is larger than the SFR(Ha) by a factor 
of 5.4 and 2.5; however, this time because the SFR(Ha) 
is underestimated. If we compare the SFR(24/im) to 
the SFR(Ha)s£;D, the SFR(24/im) is larger by factors of 
18.7, and 18.6. 

4. ARE THE BLAST-DETECTED GALAXIES TYPICAL OF 
z - 2 MASSIVE GALAXIES? 

The reasonable agreement between the FIR-derived 
and Ha derived SFRs in the previous section shows that 
provided sufficient data to robustly determine both L/ij 
and a dust corrected L(Ha), the locally-calibrated con- 
versions of Ijiii SFR and L(Ha) — > SFR appear to 
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Fig. 3. — Top Panel: SFR derived from the 24/xm fluxes and the 
luminosity dependent CEOl templates vs. the SFR determined 
from the dust corrected L(Ha) for the fuU sample of K-selected 
galaxies in Kriek et al. (2008). The luminosity dependent tem- 
plates systematically overestimate the SFR(Ha) for the majority 
of galaxies. Bottom Panel: Log-average of the SFR derived from 
the 24/im fluxes and all 105 CEOl templates vs. SFR(Ho). An 
average of the local templates provides IR SFRs that are in better 
agreement with the SFR(Ha). The luminosity dependent method 
requires that the majority of galaxies use the most luminous tem- 
plates, whereas the template averaging forces the L/^ to come 
from cooler templates. The better agreement from cooler tem- 
plates suggests the entire sample of massive z ~ 2 galaxies would 
be consistent with being "scaled up" versions of lower luminosity 
galaxies. 

provide consistent results for these galaxies. Before we 
use these two galaxies to draw conclusions about the star 
formation properties of the larger population of z 2 
massive galaxies, we must consider if they are represen- 
tative. With 24/im flux densities of 570 ± 18 /iJy and 
1417 ± 28 /iJy, these two galaxies are the brightest 24/im 
sources in the entire Kriek et al. (2008) sample. That 
they are the only ones detected in both the BLAST and 
APEX observations of the ECDFS is potentially a strong 
selection effect. 

One way to test how representative the BLAST- 
detected galaxies are of the entire K-selected sample is 
to compare the SFR(24/im) and SFR(Ha) of the full 
sample and see if it shows the same disagreement as 
the BLAST-detected galaxies. MIPS 24/im observations 
of the MUSYC flelds were obtained as part of Spitzer 
program GO-30873 (PI: Labbe). We reduced that data 
with the standard MOPEX procedures and derived 24/im 
fluxes for galaxies in the Kriek et al. (2008) sample us- 
ing the method described in § 2. There are 24 galaxies 
in the Kriek et al. sample where Ha is observable in the 
NIR window from the ground, and of these galaxies 15 
are detected at 24/xm. The MUSYC 24/im data is much 
shallower than the FIDEL data, so the non-detection of 
some galaxies, depending on their SFR and redshift is 
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expected. For the undetected galaxies we calculate 5a 
upper limits based on empty aperture fluxes in the 24/im 
mosaics. 

In the top panel of Figure 3, we plot the SFR(24/im) 
vs. the SFR(Hq;) for these galaxies. Galaxies that are 
candidates for hosting an AGN from either an X-ray de- 
tection, emission-line ratios (Kriek et al. 2007) or for 
having a powerlaw SED in the IRAC bands (Muzzin 
et al. 2009) are plotted as red symbols (this includes 
both BL AST-detected galaxies, see § 5.3). For refer- 
ence, the BLAST-detcctcd galaxies are plotted as stars 
at their SED- fit SFR(L/7j) locations with arrows con- 
necting them to their SFR(24^m). Just like the BLAST 
sources, the majority of galaxies have a SFR(24/^m) that 
is systematically larger than the SFR(Ha). This suggests 
that the local luminosity dependent templates probably 
overestimate the IR SFR of the majority of 2 ~ 2 mas- 
sive galaxies. 

Franx ct al. (2008) and Wuyts et al. (2009) also no- 
ticed that the SFRs of 2 ~ 2 galaxies tended to be larger 
when inferred from 24//,m flux and the local luminosity- 
dependent templates than the SFR determined using a 
fit to the stellar SED (SFR(SED)). They advocated using 
the L//f taken from the log-average of the local templates 
and argued that this produces SFRs that are in better- 
agreement with the SFR(SED). For comparison, we de- 
termine the SFR(24/im) using that method (although we 
use the CEOl templates, whereas they use the Dale & 
Helou 2002 templates), and compare with the SFR(Ha) 
in the bottom panel of Fig 3. 

Clearly the agreement between SFRs is much better 
when SFR(24/xm) is determined using template averag- 
ing method. Given that when using the luminosity de- 
pendent templates, most of our galaxies require the most 
luminous templates to match the 24/im flux, it suggests 
the better agreement between SFRs with the template 
averaging method is because it includes numerous lower 
luminosity templates in the average. This suggests that 
the larger sample of K-selected galaxies, like the BLAST- 
detected galaxies, can be described as scaled up versions 
of lower luminosity galaxies. Furthermore, it appears 
that despite their exceptionally large FIR luminosities, 
the cooler SEDs of the BLAST-detected galaxies may be 
representative of the larger population of massive galax- 
ies at ^; ~ 2. 

5. DISCUSSION 

5.1. Dust Temperatures and Evidence for "Scaled Up" 
Cool Galaxies 

It is remarkable that despite an Ljji ~ 10^^ Lq, both of 

the BLAST-detected z ~ 2 galaxies are bettcT fit with lo- 
cal templates that have PAH features and a cool "cold" 
dust bump. In the nearby imiverse, galaxies with L/^ 
^ 10^'^ L0 tend to have both the warmest "cold" dust 
(TdMst ~ 10 - 70K, emission peaking in the FIR) and 
significant amounts of "hot" dust (T^ust ~ 70 - 500K, 
emission peaking in the MIR). Although both of these 
dust components can come from deeply embedded star 
formation (e.g., Tran et al. 2001), it is thought that 
the dominant radiation source in most, if not all nearby 
galaxies with L ~ 10^"^ Lq is an AGN (e.g., Genzel et al. 
1998; Lutz et al. 1998). 
By contrast, the best fit templates of the 2; ~ 2 galaxies 



are much cooler ULIRG and LIRG templates which are 
most frequently associated with star forming galaxies in 
the local universe (e.g., Armus 2009). These templates 
also tend to have stronger PAH features that are more 
prevalent because the "hot" dust component from the 
AGN no longer dominates the MIR emission. 

The first claims that IR luminous galaxies at z ~ 2 may 
be scaled up versions of lower luminosity cooler galaxies 
came from 24/xm observations of submillimeter-selected 
galaxies (SMGs). Pope et al. (2006) showed that the 
SEDs of SMGs peak at longer wavelengths than expected 
based on an extrapolation of their 24/im flux using local 
templates, arguing that SMGs are cooler at a given lumi- 
nosity than the local templates. Subsequent MIR spec- 
troscopic observations of SMGs from Spitzer by Lutz et 
al. (2005), Pope et al. (2008), and Menendez-Delmestre 
et al. (2009) all found strong PAH features and relatively 
weak MIR continuum in the SMGs suggesting much less 
" hot" dust than is usually found in local galaxies of com- 
parable L/fl. Similarly, Papovich et al. (2009) found a 
weaker than expected rest-frame 24/xm luminosity for a 
lensed SMG, arguing for the lack of a "hot" dust compo- 
nent, and hence a scaled up cool galaxy. 

With access to the FIR we can determine a dust tem- 
perature (T(j) independent of the MIR SED and com- 
pare this to the local Td-Ljji relation. We fit the 160/im 
— > 870/im data (50/im — > 300/um rest-frame) to modi- 
fied blackbody curves of the form S(j^) = Aiy'^B{iy,Td), 
where A is a normalization, B is the Planck function, 
and accounts for frequency-dependent dust emissiv- 
ity. Even with detections in numerous bands the S/N of 
the data is too low to constrain both j3 and T^ simul- 
taneously so we have assumed /3 = 1.5, a typical value 
for submillimeter galaxies. The 70/im points are omit- 
ted in the fitting because they correspond to rest-frame 
23/im and are too far into the MIR to include in a single- 
component blackbody fit. We note that including the 
70/im data actually produces identical best-fit temper- 
atures; however, the of fits are substantially larger 
because no single-temperature model can describe the 
23/im — 300/im rest-frame simultaneously. 

We find dust temperatures of 40l^K and 41^7K for 
ECDFS-4511 and ECDFS-12514, respectively. These 
correspond to rest-frame monochromatic IRAS color 
measurements, C = logio(f6o/im/fioo/im), of -O.OSto og 
and -0.03^0.20' respectively. Comparing these to the C- 
L/_R relation measured by Chapin et al. (2009) in the 
nearby universe using IRAS shows that these galaxies 
have temperatures typical of nearby galaxies with Lm 
= 10" - 10^2 Lq. They are also similar to the T^ 
of local LIRGS/ULIRGS measured by Clements et al. 
(2010). Therefore, with the combined Spitzer, BLAST 
and APEX data we can now confirm that some massive 
galaxies at z ^ 2 not only have PAH-dominated MIR 
SEDs similar to lower luminosity local galaxies, but they 
also have "cold" dust SEDs consistent with being scaled 
up versions of lower luminosity galaxies. This suggests 
that star formation may be occurring in a different envi- 
ronment within galaxies at z ~ 2 compared to z ~ 0. 

Local galaxies with SFRs > 100 Mq yr~^ tend to 
have compact nuclear starbursts which are the caiise of 
the warmer dust temperatures (e.g., Sanders & Mirabel 
1996). The z ~ 2 galaxies have SFRs ~ 1000 Mq yr"!, 
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but similar dust tcmpcratiircs which suggests that the 
star formation may be occurring in multiple large star- 
burst regions throughout the galaxy, rather than concen- 
trated in a single compact starburst which would pre- 
sumably have a substantially hotter dust temperature. 
Indeed, a Icnsed z > 2 galaxy with precisely such char- 
acteristics has recently been observed by Swinbank et 
al. (2010a) using interferometric submillimeter observa- 
tions. 

Kriek et al. (2009a) presented SINFONI-IFU Ha spec- 
tra of both these galaxies. Their maps show that the Ha 
emission is resolved on scales of a few kpc, and is also 
consistent with the extended star formation hypothesis. 
The greater availability of gas at high redshift may per- 
mit more spatially extended star formation at a high rate, 
something that is uncommon in local galaxies. 

5.2. Morphologies 

In the inset of Figure 1 we show 3" x3" NIC2 F160W 
images of the galaxies from Kriek et al. (2009a). At 
z 2 these correspond to rest-frame V-band morpholo- 
gies and span an angular size of ~ 25 kpc. Like most 
of the high redshift ULIRGs and HLIRGs previously ob- 
served with HST (e.g., Chapman et al. 2003; Dasyra et 
al. 2008; Swinbank et al. 2010b) the galaxies clearly 
have disturbed morphologies. Kriek et al. (2009a) mea- 
sure effective radii (Re) of 5 and 3 kpc for ECDFS-4511 
and ECDFS-12514, respectively, which means that both 
galaxies are comparable in size to local galaxies of simi- 
lar mass. Rest-frame optical morphologies make it diffi- 
cult to identify the location of the star formation (which 
is presumably highly obscured); however, the disturbed 
nature of both galaxies shows that these are complex sys- 
tems, possibly "stream- fed" galaxies (Dekel et al. 2009) 
or mergers. Unlike the FIR SEDs, the morphologies do 
not appear to be scaled up versions of local star forming 
disk galaxies. 

5.3. AGN Contribution to the MIR and FIR Emission 

Despite SEDs that are associated with local star form- 
ing galaxies, part of the MIR and FIR emission in both 
galaxies almost certainly comes from an AGN. Using 
emission line ratios and the BPT diagram (Baldwin et 
al. 1981), Kriek et al. (2007) found that both galax- 
ies were in the region of the diagram that Kewley et al. 
(2006) considered to have "composite" spectra with con- 
tributions from both star formation and an AGN. 

ECDFS-4511 is detected in the 2 Ms Chandra map of 
the CDFS (Luo et al. 2008), and has an ~ 2x10''"' 
erg s~^, implying it hosts a Seyfert-luminosity AGN. 
ECDFS-12514 is located in the 250 ks map of the ECDFS 
(Virani et al. 2006), but is imdctccted. Although unde- 
tected in the X-rays, ECDFS-12514 is formally classified 
as an IRAC powerlaw galaxy (F^ cx ly", a < -0.5; e.g., 
Donley et al. 2007, see Figure 1) which makes it a can- 
didate for an obscured AGN. 

Without more data we cannot constrain the precise 
fraction of the MIR and FIR light that comes from the 
AGN; however, the consistency of the MIR/FIR SED 
with local templates that are star formation dominated, 
and the "composite" emission line ratios suggest that 
young stars are probably the dominant energy source 
heating the dust. Pope et al. (2008) and Menendez- 
Delmestre et al. (2009) showed that the MIR SEDs of 



most SMGs (the BLAST-detected galaxies are formally 
SMGs) are still dominated by PAHs and that AGN only 
contribute ^ 30% of the total MIR emission. Putting the 
evidence together the most plausible scenario is that both 
galaxies host a moderate-luminosity AGN, but the AGN 
does not dominate the MIR/FIR SED. Nevertheless, we 
note that both our SFR(L/ij) and SFR(Ha) should be 
formally considered upper limits as the AGN will con- 
tribute to both of these. 

5.4. Implications for SFRs at z ^ 2 and MIR-excess 

Galaxies 

The results of our FIR SED fitting show that caution 

is required when determining and interpreting the SFRs 
of z ^ 2 galaxies based on limited data. Using only 
24/im photometry to infer SFRs at z ^ 2 may work well 
for lower luminosity, or less obscured galaxies (e.g., Pa- 
povich et al. 2007; Reddy et al. 2010); however, the 
local templates can overestimate the SFR(L/i{) of the 
most FIR luminous galaxies by as much as a factor of ~ 
6 (sec also Papovich et al. 2007; Murphy et al. 2009). 
Likewise, the Balmer decrement measurements for these 
galaxies seem to indicate extra extinction toward the star 
forming regions. It appears that on average this relation 
is similar to that suggested by Calzetti et al. (1997); 
however, the lack of secure Hfl detections in our sample 
precludes quantitative statements. To correctly deter- 
mine the SFR(Ha) of individual galaxies, rather than 
ensemble averages requires high quality Balmer decre- 
ment measurements. 

The extra extinction towards the star forming regions 
suggested by the Balmer decrements also indicates that 
the UV SFRs for some of these IR luminous galaxies 
may not be reliable because the A^, is large enough that 
parts of the star forming regions are optically-thick. The 
SFRs estimated from a fit to the stellar populations 
(SFR(SED); effectively a dust corrected UV SFR) by 
Muzzin et al. (2009) for ECDFS-4511 and ECDFS-12514 
are 168163 yr^^ and U5tl^° ^0 Y^'^' and under- 
estimate the SFR(L/fl:) of the galaxies by factors of 5.2 
and 3.5, respectively. If we compare the SFR(SED) to 
the SFR(24^m) as in e.g., Daddi et al. (2007b) we find 
that these galaxies have MIR SFRs 18 25 times greater 
than the UV SFRs. and are extremely strong "MIR ex- 
cess" galaxies. 

The purpose of the Daddi et al. (2007b) classifica- 
tion was to identify candidate obscured AGN, and given 
that it is likely that both galaxies contain an AGN it 
appears that method successfully identifies such sources. 
However. MIR excesses of factors of 18 25 would sug- 
gest that the AGN is the dominant MIR energy source 
for both galaxies, something that is unlikely given the 
shape of the MIR/FIR SEDs and the emission line ratios. 
This calls into question whether galaxies with weaker 
inferred MIR excesses, such those near the cutoff of 
SFR(UVcorr)/SFR(MIR-hUV) ~ 3 suggested by Daddi 
et al. (2007b) have a significant contribution from an 
AGN. 

5.5. Comparison to Models 

Our results suggest that a substantial population of 
massive z ~ 2 galaxies are "scaled up" versions of lower 
luminosity galaxies. Clearly such a population needs to 
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be understood within the context of modern galaxy for- 
mation and evolution models. Recently, Hopkins et al. 
(2010) have used their merger-driven evolutionary model 
to predict the IR luminosity functions of galaxies and 
quasars between < 2: < 6. In their model they ar- 
gue that there is a threshold in hjR above which disks, 
merger-driven star formation, and obscured AGN are the 
dominant IR energy source of galaxies. At 2: = these 
limits are log(L/fl;) > 11.5 for merger-driven star forma- 
tion, and log(L/i^) > 12.5 for obscured AGN. Based on 
the evolution of the merger rate and gas fraction to z ^ 
2 they suggest that these thresholds evolve in log(L/i^) 
to log(L7/j) > 12.5 for merger-driven star formation and 
log(L7jj) > 13.0 for obscured AGN. With log(L/H) ~ 
13.0, the BLAST galaxies lie at an L/jj that is near the 
suggested transition point between merger-driven star 
formation dominated or AGN dominated. 

Interestingly, the FIR SEDs of these galaxic^s do not 
look like those of the AGN-dominated HLIRGs detected 
by ISO (Genzel et al. 1998), or of the average nearby 
ULIRG, which are mostly merger-driven starbursts (e.g., 
Sanders et al. 1998). Instead, these galaxies look like 
scaled up versions of even lower luminosity systems, 
which are not universally associated with mergers or 
AGN. Taken at face value, it appears they do not fit 
well within the Hopkins et al. (2010) description. 

On the other hand, the NIC2 morphologies of the 
galaxies are clearly disturbed enough to be considered 
candidates for undergoing mergers (Kriek et al. 2009b). 
If these galaxies are merger-driven starbursts, it may be 
that they are missing the centrally-concentrated star- 
burst which produces the much hotter dust seen in these 
galaxies in the local universe. As Hopkins et al. (2010) 
point out, the efficiency of a merger at funneling gas and 
stars into the core of the galaxy is proportional to the 
gas fraction, with higher gas fraction mergers being less 
efficient. If both galaxies have high gas fractions most 
of the star formation could be occurring on much larger 
physical scales than similar mass starbursts in the local 
universe, hence the cool dust SED. 

Galaxies with similar stellar masses (M ~ 10^^ M©) 
and gas fractions as high as 80% have been observed in 
the distant universe (e.g., Tacconi et al. 2010) and could 
be analogues of the BLAST-detected galaxies. If that 
is the case, it begs the question of what are the descen- 
dants of extended, massive, gas-rich star forming galaxies 
at z ^ 21 Local spheroids seem like an obvious candi- 
date; however, some fraction of the progenitors of those 
systems are clearly small and quiescent at z ~ 2 (e.g., 
van Dokkum et al. 2010; Hopkins et al. 2010; Bezanson 
et al. 2009; van der Wei et al. 2008). If the large and ac- 
tive BLAST-sources are also progenitors of massive local 
spheroids it suggests there may be very different evolu- 
tionary paths which arrive at the same galaxies. 

These galaxies also appear to be consistent with the 
"stream-fed" galaxies proposed by Dekcl et al. (2009) 
and the scaled up "turbulent disks" discussed by Genzel 
et al. (2008). In their model, Dekel et al. (2009) suggest 
that massive galaxies may grow quickly through rapid 
cold gas accretion and clumpy star formation which fun- 
nels stars into a central bulge. Given that these cold flows 
often separate into multiple large star forming clumps, it 
might be expected that these clumps produce cooler dust 



temperatures at fixed SFR than centrally concentrated 

starbursts. 

The SFRs of the BLAST galaxies are higher than the 
100 - 200 M0 yr~^ predicted for "stream- fed" galaxies 
by Dekel et al. (2009). The BLAST galaxies are the 
most IR luminous of the K-selected sample of Kriek et 
al. (2008), so if they are "stream- fed" galaxies, they may 
just represent the most active tail of the distribution. 

High quality kinematic data could potentially distin- 
guish between merger-dominated or accretion-dominated 
systems. These galaxies do have SINFONI IFU observa- 
tions; however, these data lack the S/N to reliably deter- 
mine V/(J. Therefore, at present based on the MIR/FIR 
properties and morphologies we cannot distinguish be- 
tween evolutionary models such suggested by both Hop- 
kins et al. (2010) and Dekel et al. (2009). 

6. SUMMARY 

Using a combination of data from the Spitzer, BLAST 
and APEX observatories we have presented some of the 
first well-sampled FIR SEDs for two z ^ 2 galaxies with 
rest-frame optical spectroscopy. In both galaxies we see 
a clear peak in the FIR emission at ~ 300/xm which 
is caused by the redshifted thermal emission from dust. 
The FIR SEDs of the galaxies are well described by lo- 
cal templates; however, the L/^ appears to be an order 
of magnitude larger than the associated with the 
locally calibrated templates. This suggests that these 
galaxies are likely to be scaled up versions of "cool" local 
galaxies, rather than analogues of nearby HLIRG tem- 
plates. 

A comparison of SFR(L/i^) and SFR(Ha) for these two 
galaxies shows that they are consistent; however, if we 
use only the 24/im photometry and the local luminosity 
dependent templates to determine the SFR we find that 
this overestimates the SFR(Ha) by as much as a fac- 
tor of 6. A larger sample of 24 galaxies from the Kriek 
et al. GNIRS sample shows the same relation between 
SFR(24^m) and SFR(Ha) and suggests that many mas- 
sive galaxies at z ^ 2 arc likely to be scaled up versions 
of lower luminosity galaxies. 

Although based on just a handful of galaxies, our re- 
sults show that FIR luminous galaxies at z ~ 2 are clearly 
different than their counterparts in the local universe. 
In the coming years the Herschel and JWST observa- 
tories will provide us with large amounts of high-quality 
FIR observations on the population of luminous galaxies 
at z >~ 2. Likewise, upcoming NIR multi-object spec- 
trograph will provide better Ha and H^ measurements 
for those galaxies, and the ALMA observatory will pro- 
vide spatially-resolved FIR maps of these galaxies. This 
plethora of new, high-quality data will allow us to bet- 
ter understand how the population of luminous, "cool" 
star forming galaxies fit into the overall picture of galaxy 
evolution. 

We thank the anonymous referee whose comments 
helped improve the quality of the paper. Support for pro- 
gram HST-GO-11135.08 was provided by NASA through 
a grant from the Space Telescope Science Institute. 
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